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Individual processive molecular motors, of which conventional kinesin is the most studied quan- 
titatively, move along polar molecular tracks and, by exerting a force F = (F x ,F y ,F z ) on a tether, 
drag cellular cargoes, in vivo, or spherical beads, in vitro, taking up to hundreds of nanometer-scale 
steps. From observations of velocities and the dispersion of displacements with time, under mea- 
sured forces and controlled fuel supply (typically ATP), one may hope to obtain insight into the 
molecular motions undergone in the individual steps. In the simplest situation, the load force F may 
be regarded as a scalar resisting force, F x < 0, acting parallel to the track: however, experiments, 
originally by Gittes et al. (1996), have imposed perpendicular (or vertical) loads, F z > 0, while 
more recently Block and coworkers (2002, 2003) and Carter and Cross (2005) have studied assisting 
(or reverse) loads, F x > 0, and also sideways (or transverse) loads F y 7^ 0. 

We extend previous mechanochemical kinetic models by explicitly implementing a free-energy 
landscape picture in order to allow for the full vectorial nature of the force F transmitted by the 
tether. The load-dependence of the various forward and reverse transition rates is embodied in 
load distribution vectors, 6^ and 0~ , which relate to substeps of the motor, and in next order, in 
compliance matrices 77+ and rjj . The approach is applied specifically to discuss the experiments of 
Howard and coworkers (1996) in which the buckling of partially clamped microtubules was measured 
under the action of bound kinesin molecules which induced determined perpendicular loads. But in 
the normal single-bead assay it also proves imperative to allow for F z > 0: the appropriate analysis 
for kinesin, suggesting that the motor "crouches" on binding ATP prior to stepping, is sketched. It 
yields an expression for the velocity, V(F X , F z ; [ATP]), needed to address the buckling experiments. 



INTRODUCTION 



A processive motor protein Q] is an individual protein 
molecule that in an appropriate aqueous solution binds to 
a linear periodic directed molecular track and, when fu- 
eled via diffusion by suitable molecules (specifically ATP 
in cases of most interest), takes tens to hundreds of dis- 
crete steps along the track before dissociating. Such a 
motor steps overwhelmingly in a single characteristic di- 
rection, which we will identify as parallel to the positive 
x-axis: see Figure^ Thus conventional kinesin walks to- 
wards the plus or fast-growing end of a microtubule while 
myosin V moves towards the plus or 'barbed' end of an 
actin filament 1]. In stepping along its track a motor 
exerts a tension F — (F x ,F y ,F z ) on a molecular tether 
the other end of which is, in vivo, bound to some cellu- 
lar organelle or vesicle while for in vitro experiments it 
is firmly attached to a silica or plastic bead that may be 
controlled with the aid of an optical trap For most 
purposes one may regard the molecular track as fixed in 
space: relatively rigid microtubules are typically clamped 
to a microscope slide while the more flexible actin fila- 
ments can be stretched between two further beads that 
are attached to the filament ends and held in place by 
individual optical traps . 

In the conceptually most straightforward experiments 
HHUm a motor (M in Fig.^l advances along the track in 
a positive direction in a solution of fixed fuel concentra- 
tion, say [ATP], under a resistive load the ^-component 
of which, F x < 0, is observed or controlled. Measure- 
ments of the mean velocity and the dispersion 0, Q : 

V « (x(t))/t, and D a [(x 2 (t)) - (x(t)f]/t, (1) 




FIG. 1: Schematic representation of a motor protein, M, 
bound to a molecular track (depicted as a microtubule with + 
and — ends) . The point of attachment of the tether to the mo- 
tor body is marked P. The force F = (F x , F y ,F z ) denotes the 
tension transmitted by the tether. The cartesian coordinate 
system has its origin at the point of attachment of a station- 
ary motor at a binding site further along the track. The sense 
of the mean velocity, V, under low loads is indicated by an 
arrow parallel to the track. 



as functions of F x and [ATP], where x(t) represents the 
displacement of the motor along the track at time t, may 
then be analyzed 0, 0, with the aim of extracting 
details of the mechanism by which the motor takes indi- 
vidual steps. In particular one would like to identify and 
quantify any substeps or intermediate motions. 

To investigate further experimentally it is rather nat- 
ural to impose an oppositely directed reverse or assisting 
load on the motor so that F x > 0. Indeed, such experi- 
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ments have been performed on kinesin, first by Coppin et 
al. 0, and, more recently, by Block and coworkers [3. Ho| . 
and by Carter and Cross In a preliminary compu- 
tation to gain theoretical insight regarding the effects of 
reverse loading, Fisher and Kolomeisky (FK) exam- 
ined the fairly noisy data of Coppin et al. ||| for the mean 
velocity of kinesin under assisting loads up to F x = + 6 
pN. Merely by continuing the formulae developed for neg- 
ative, i.e. resistive values of F x , analytically through 
F x = to positive values, FK obtained an apparently 
reasonable description of the observed 'acceleration ra- 
tios,' V(F X = 5 pN)/V(F x = 0). These varied from about 
3.0 to 1.4 as [ATP] increased from 5 /iM to 1 mM. The 
more recent experiments, however, do not confirm these 
results [1113, [ill; on the contrary, for [ATP] > 50 (jM 
the data indicate no acceleration or even a deceleration, 
for F x up to 10-15 pN, i.e., V(F X > 0) < V(F x =0). 

But, irrespective of the experimental observations, 
an examination of the geometry entailed in a standard 
single-bead assay with an assisting load, as illustrated in 
Figure , reveals that switching from a resistive to 



To some degree, however, the issue can be addressed ex- 
perimentally by employing beads of varying diameters. 
Thus, by reference to Figure [3 one sees that increasing 
the bead radius, R, increases F z relative to a controlled 
value of F x via 



— = cot e = 

F x 



R + Az 



[(I - Az)(2R + l + Az)} 1 / 2 



in which I is the length of the tether while Az 2 represents 
the mean squared thermal fluctuation of the bead in the 
z direction |12j . 

A decade ago, however, Howard stressed the vecto- 
rial nature of the load and the significance of indepen- 
dently observing the perpendicular or vertical compo- 
nent, Fx = F z , and measuring its effects on V. To this 
end Gittes, Meyhofer, Baek and Howard (GMBH) [b| 
devised an ingenious microtubule (MT) buckling experi- 
ment. In their protocol the minus end of an MT is chem- 
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FIG. 3: Arrangement of the microtubule (MT) buckling ex- 
periment [13] showing an originally straight MT (dashed lines) 
buckled by a kinesin motor that bound to the MT and pro- 
ceeded to step towards the plus end at a velocity V that de- 
pends on the varying force F transmitted by the tether. 



FIG. 2: Depiction of the configuration of a bead, of radius 
R, and a tether, of length I, attached to the motor body at 
P when the motor (moving on a microtubule clamped to a 
substrate) is subjected to an assisting, F x > 0, load. The 
offset Az denotes the root mean square thermal fluctuation 
of the bead as limited by collision with the track surface [12]. 

an assisting load should not be described merely by the 
change in sign of a scalar load force. Rather the true 
vectorial character of the force F, acting at the point of 
attachment, P, of the tether to the body of the motor 
should be recognized. Even if sideways, F y ^ 0, com- 
ponents of F may be neglected — although Block and 
coworkers 0, E3] have imposed transverse loads — a sat- 
isfactory description of the motor operation should seek 
to understand the mean velocity V, and likewise the dis- 
persion D, as functions of F x and F z . 

Clearly, it would be advantageous experimentally to 
vary F z independently of F x . With current optical trap 
technology this is not readily accomplished (since the 
traps are rather soft along the perpendicular or z-axis). 



ically attached to a glass substrate while the plus end of 
the MT can interact freely with a sparse field of kinesin 
motors tethered randomly to the same substrate. Then, 
as the free end of the MT diffusively wanders across the 
surface, a single kinesin molecule will occasionally en- 
counter the MT and bind to it. In the presence of ATP 
it then proceeds to step along the MT towards the plus 
end: see Fig. As the motor moves, it will eventu- 
ally start to buckle the MT. By recording the buckling 
process and analyzing the shapes of the MT, the two 
components of the force, F — (F X ,F Z ), in the plane of 
the substrate can be determined. Likewise the varying 
speed, V(t), along the MT can be measured. The exper- 
iment is not easy and the resulting data are quite noisy: 
nevertheless, GMBH felt able to conclude that the ver- 
tical component, F z , resulted in an increase in the overall 
mean speed. 

In this article we extend the previous quantitative 
mechanochemical kinetic approach 0, 0, Q in order to 
explicitly account for the vectorial character of the load 
force in analyzing the stochastic motions of motor pro- 
teins. The resulting formalism has been applied to the 
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recent experiments of Block and coworkers on kinesin 
[Til PHI . Contrary to their model description, which im- 
plies sideways (or transverse) motions of the motor while 
stepping, we find no cause to invoke displacements of the 
point of attachment outside the (x, z) plane. However, 
our analysis indicates that a kinesin motor 'crouches' on 
binding ATP, that is, the point of attachment of the 
tether moves downwards towards the microtubule track, 
apparently by about 0.5 to 0.8 nm before rapidly com- 
pleting a unitary 'swing' or diffusive step of close to 8.0 
nm transferring the motor to the next binding site on the 
track [HD3. 

Using the corresponding fitted kinetic parameters, we 
have revisited the Gittes et al. microtubule buckling ex- 
periment |l3[ . As we explain below, our analysis avoids 
a simplifying assumption made in their discussion and 
seems qualitatively consistent with the buckling data; 
however, it indicates that V(F X ,F Z ; [ATP]) generally de- 
creases when F z increases (at fixed F x and [ATP]) rather 
than increases as GMBH argued. Nevertheless, one may 
note that in the buckling experiments the motor moves 
progressively away from a curved and, hence, stressed 
region of the microtubule. This raises the interesting 
question as to the degree to which the curvature of a 
microtubule protofilament might affect the motility of 
kinesin. 



2. INTERMEDIATE STATES AND SUBSTEPS 

The modeling of motor protein motility may be carried 
out at different levels P, E3, HE OH Hi! • For concreteness 
and relative simplicity, we will focus on motors powered 
by the hydrolysis of ATP and moving like kinesin via 
steps of fixed size. Then there are specific binding sites 
located at positions x = Id (I = 0, ±1, ±2, • • • ) along the 
periodic track where d is the step distance corresponding 
to the track periodicity. For a microtubule one has d ~ 
8.2 nm, representing the size of a tubulin dimer. In the 
absence of molecular fuel, [ATP] = 0, the motor will be 
bound at a specific site / in a nucleotide-free state to be 
labeled 0;. We will overlook the retention of ADP in 
the weakly bound head of kinesin 0, 0] and neglect the 
rates of spontaneous dissociation from the track: but see 

When, accepting the evidence for "tight coupling" 
the motor translocates from site I to the adjacent site 
I + 1, by processing a single fuel molecule, it undergoes, 
in the simplest linear reaction sequence considered here, 
a series of N intermediate mechanochemical transitions 
from states j = 0; to 1; to • •• to (N — 1); to Ni = 0i+i- 
When a motor is powered by the hydrolysis of ATP, the 
biochemical evidence indicates that N — 4 distinct enzy- 
matic states should be recognized 0; however, the de- 
gree to which these are linked to significantly different 
mechanical states is a significant object of investigation. 
Clearly, the simplest nontrivial model requires N = 2 
states: the transition from state 0/ to 1; then corresponds 



to the binding of ATP which is followed by the hydroly- 
sis process that is completed, with loss of ADP and Pi, 
when the motor moves on to state 2/ = 0;+i. Again for 
simplicity, we overlook the coupling between two distinct 
heads of a motor as entailed in the so-called hand-over- 
hand motion now well established for kinesin: see, e.g. 
[2TL I2H I24L l2fij . Nevertheless, this can be accommodated 
formally simply by doubling the nominal step size d and 
allowing for twice as many intermediate states to describe 
periodic double-steps. (See, e.g., the analysis used in Q 
for myosin V which, indeed, exhibits steps of fluctuating 
size.) 

In the simplest, traditional chemical kinetic pictures 
Plr^] one introduces transition rates, forwards and back- 
wards, between states along a reaction pathway: thus we 
will associate rates Uj and Wj with the transitions from 
state ji forward to state (j + 1)/ and backward to state 
(j — 1);, respectively. Note that owing to the periodicity 
of the stepping process, the rates Uj and Wj are indepen- 
dent of I. 

However, to describe the action of motor proteins un- 
der variable loads it is essential to describe the depen- 
dence of the various rates on the loads imposed, or, 
complementarily, on the stresses induced: and how the 
influence of the load is distributed over the different 
mechanochemcial transitions should be a prime question 
[j| . Since purely chemical steps are typically very fast on 
the scale of the overall mechanical movements, a rather 
basic theoretical picture ^E] postulates a distinct free 
energy surface, depending on molecular coordinates, for 
each separate mechanically fluctuating biochemical state. 
Motions occur by diffusion through the multidimensional 
configurational space while the probability of a formally 
instantaneous "chemical jump" from one surface to an- 
other varies with the configuration. One may argue, how- 
ever p| , that for practical purposes this scheme may be 
reduced to an effective mechanochemical description in 
which points along the traditional chemical "reaction co- 
ordinate" are, for a motor protein, simply identified with 
specific positions of the motor along the linear track. 

By this route, as previously adopted [a,!?!, one may, 
indeed, hope to identify substeps, do, d\, in which 
the motor moves from a binding site, say, at xq to the 
next site xo + d via intermediate locations x\ = xq + do, 
X2 = + dg + d\ , ■ ■ • , with ^2 j—o 1 dj = ^- Likewise, 
within the chemical picture, one may identify successive, 
unstable "transition states," say j ; + = j7+i> along the 
reaction coordinate each lying between the locally stable 
intermediate states ji and In mapping the reaction 
coordinate onto the track, one may then decompose a 
substep dj from Xj to Xj+i into the sum (dj + dj +1 ) 
thereby locating the transition state at x^ = xj + d^ = 
Xj + i — dj +1 = xj +1 : see Fig. ^ of for depictions of 
such mappings with N = 2 and N = 4, as fitted to data 
for kinesin under resisting (F x < 0) loads 0, Q . 

It is rather clear that in such a treatment, in which the 
reaction pathway is assumed to be parallel to the £-axis, 
only the longitudinal component, F x , of the load will be 
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coupled to displacements of the motor and so affect the 
rate constants. However, the initial substep predicted by 
this approach (do — 1.8 -2. f nm) appears to be inconsis- 
tent with high-resolution spatio-temporal observations of 
individual and averag ed steps [H M 111 El 

Further- 
more, this approach has failed to account satisfactorily 
for velocity measurements of kinesin under reverse load- 
ing [U lid Hl| . Accordingly, apart from general theo- 
retical considerations, it seems important to extend the 
previous discussions. 

Evidently, a significant next step in considering the 
load-dependence is to recognize the vectorial character 
of the imposed force by allowing states along the reac- 
tion pathway to be mapped onto movements of the motor 
within the full three-dimensional space of the track, mo- 
tor, tether, and cargo: see Fig. ^ More concretely we 
may hope to follow the course of the point of attachment 
P => r(t) = [x(t),y(t), z(t)] of the tether to the motor as 
steps are taken: see Figs. ^ and [3 As varying values, z(f) 
and y(t), arise one can view the motor as moving up or 
down or swinging from side to side when it walks along 
the track from a binding site I to the next binding site 
I + 1. The simplest (N = 2)-state model embodying this 
concept is illustrated in Fig. 0] It has been supposed, 




FIG. 4: Schematic depiction within an (N = 2)-state model 
of the motion in the (a;, z) plane of the point of attachment, 
P = r = (x, y, z), of the tether to the motor body. Two bound 
nucleotide-free states, labeled and = 2, are shown to- 
gether with a single (N — 1 = 1) intermediate mechanochem- 
ical state, 1, and the two (JV = 2) corresponding transition 
states (0 + = 1~ and l + =2~). See further discussion in the 
text. 

first, that the motion of P may be regarded as confined 
to the (%, z) plane, i.e., that sideways (or transverse) mo- 
tions may be neglected. Of course this should be subject 
to experimental test: but for kinesin [UEJEI it proves 
quite adequate |l5j . 

Then in Fig. 0] the rectangular boxes labeled "0" and 
"2 = 0" represent bound, nucleotide-free states at suc- 
cessive sites I and I + 1, along the track. The circle la- 
beled "1" denotes the ATP bound state; this state may 
also include subsequent ADP-Pi and ADP biochemical 



states. This TV = 2 level model, however, cannot dis- 
tinguish situations in which the hydrolysis, ADP, and 
Pi-release chemical reactions take place physically at (or 
close to) the mechanical state "1" or close to "2" . (And 
recall also the remark above concerning hand- over- hand 
motion.) 

The (unlabeled) crosses in Fig.0|represent the location 
of the two transition states, + = l - and 1 + = 2~. 
Also introduced in the figure are the dimensionless load 
distribution factors 8j + , , 8j + and 6^~~ (following 
5, 6]). These turn out to be crucial fitting parameters 
in describing experimental data. Evidently they serve to 
specify the substeps via 



9fc 



and 



dj = 



+ d-+!)d, (3) 



while vertical (or perpendicular) displacements of the 
point of attachment are specified by 



Azf = efd 



and 



Az j = (q + + ^ 1 )d, (4) 



with, of course, similar expression for Ay^ in terms of 
#V . Since the stepping is periodic one must finally have 



JV-l 



3=0 



N—l 

£( 

3=0 



9f 



) = o. 



(•5) 



LOAD-DEPENDENCE OF TRANSITION 
RATES 



Now, under any external load, F, the rate constants, 
Uj and Wj, must change. But how? To answer, let us 
accpet a set of "free-energy landscapes" , for distinct bio- 
chemical states each depending on some set of molecular 
configurations 18]. Then one can contemplate a "re- 
duced landscape" in which the set of distinct landscapes 
is, in effect, combined by identifying a traditional reac- 
tion coordinate, while at the same time retaining a par- 
ticular mechanical coordinate, like the displacement, x, 
of the motor along its track 20]. (One may, incidently, 
note a critique |2(j of this general approach.) It is nat- 
ural in the present case, as we have demonstrated, to 
extend the retained mechanical coordinate to be the vec- 
tor r specifying, in real space, the point of attachment, 
P, of the tether to the motor. Then, if one supposes the 
reaction pathway can, at least on average, be mapped on 
to the motion of P, one may dispense with any explicit 
consideration of the reaction coordinate. Thus we are 
led to postulate a more-or-less well defined, presumably 
smooth and diffcrcntiablc free-energy function, say <I>(r), 
which in the absence of a load (i.e., F = 0) respects the 
periodicity of the track so that 



$(r) = <f>(r + dx), 



(6) 



where x is a unit vector parallel to the x axis. (See also 
Chaps 15, 16, etc., and Sachs and Lecar [27j].) 
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The various mechanochemical states 0/, 1;, ■ ■ ■ , ji, 
are then to be identified with corresponding valleys or po- 
tential wells, i.e., minima in $(r), located, say at rj+ldx 
for j = 0, 1, • • • (mod N) and I = 0, ±1, ±2, • • • . Follow- 
ing the traditional chemical picture, successive valleys 
along the reaction path will be linked via free energy 
barriers corresponding to the respective transition states. 
These, in turn, will be represented, by cols (or passes or 
saddles) in <I>(r) at points r+ = rj +1 . 

When a force F is exerted on the motor's tether the 
(now reduced) free energy landscape must be deformed. 
It is then most reasonable (but surely not fully 'realistic') 
to suppose that 

$(r;F) = $(r) - F-r. (7) 

Now, once again in accord with traditional chemical 
reaction rate theories (see, e.g. 0), the rate Uj of the 
forward reaction from a state j to j + 1 will be dominated 
by the Boltzmann factor for overcoming the correspond- 
ing barrier; and the same goes for the backwards reaction 
to j — 1. Thus, if <&j(F) is the free energy level at the 
bottom of the jth valley while $>~j(F) and $j(F) are 
the height of the associated forward and rearward col, 
respectively, we may suppose 

Uj (F) oc e-l*tW-*iWW* T , 

Wj (F) oc e -[*7W-*^)]/ fc BT (g) 

To proceed so as to obtain the force-dependence of the 
rates up to quadratic order in F on the basis of this 
landscape picture, consider first the forward rate con- 
stant Uj. In the absence of the load the corresponding 
valley is located at rj while the appropriate col (describ- 
ing the transition state j + ) is at r+ = rj + Oj~d where 
we have introduced the forward load distribution vector 
0+ = (#J + , Expansion of the free energy in 

the valley then yields 

*(r) = ^ + i(r - r^-A^r - r,) + 0(\r - r,\ 3 ), (9) 

where $j = &(rj) while Aj is a positive definite sym- 
metric 3x3 matrix with elements A-^ = (d 2 ^/dXdfi) 
evaluated at r,-, where A, /i = x, y, z. On the other hand, 
one may expand the free energy around the col in the 
form 

$(r) = $t + i( r - r+)-A+-(r - r+) + ■ ■ • , (10) 

where = &(r~j~) and A+ is again a 3 x 3 matrix with 
elements (<9 2 <!>/<9A<9/i) but now evaluated at r^. Note 
that by the character of a col or saddle point, the ma- 
trix A* 1 ", which is identical to AJ +1 , must be indefinite 
with at least one negative eigenvalue. The correspond- 
ing eigenvector identifies the optimal direction of the re- 
action, projected into r-space as the system enters and 
leaves the transition state. It will be appropriate for us 
to assume that the remaining eigenvalues are positive. 



The rate constant m° for F = is then proportional to 
exp(-A$+°/fc B T) where A$+° = $(r+) - $(rj) is the 
barrier height. Under a vectorial load F, however, the 
positions of both valley and col shift, the changes being 
proportional to F in leading order. The new positions 
may be found by solving the equations V$(r;.F) = 
using (JjJJ, (|10|l and J7J. Thus one finds the free energy 
difference between the shifted valley and col to be 

A$+(F) = $f(F) — $j(F), 

= A$+° - F-{r+ - rj ) - ^F-q+-F + 0(F 3 ), 

(11) 

in which appears the matrix 

Vt = (A+)- 1 - (A 3 )-\ (12) 

where by our construction, the inverse matrices of Aj and 
At are well defined. Note, again, that r+ — rj = O^d. 

By the same arguments, the free energy barrier for the 
reverse rate Wj is determined by the transition state j~ 
located at rj = rj — Ojd and is given by 

A$t(F) - $J(F)-$j(F), 

= A$T° + F-6jd+ \F-r)--F + 0(F 3 ), 

(13) 

where dj = (0j~ , ffV~ , 6j~) is the reverse load distribu- 
tion vector while 

rij = (A/)- 1 - (A,-)" 1 , (14) 

in which AJ = A^_ 1 is the matrix of the second deriva- 
tives of $ evaluated at rj = r^_ 1 . 

Finally, we may express the load-dependence of the 
rate constants as 

Uj (F) = u]e^>{ + [e+-Fd+\F-r,+ -F + 0{F 3 )]/k B T}, 

(15) 

Wj (F) = w Q 3 exp{-[0J-Fd+ \F-qJ-F + 0(F 3 )}/k B T}. 

(16) 

Evidently, the load distribution vectors 9j and 8~ sim- 
ply generalize the previous load distribution scalars 0, IE 
|7j and serve to locate the intermediate mechanochemical 
states and determine the vectorial character of the sub- 
steps linking them. By the underlying periodicity they 
must satisfy 

AT-l 

J2(et + oj) = x, (17) 

3=0 

which simply summarizes pj and J^J. 

Similarly, the matrices tj^" and rfj serve to measure 
the relative compliances of the various transition states. 
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We may note that, by (|12|) and (|14f> . the consecutive 
differences 

<-r, i+1 = (A, +1 )- 1 -(A,)- 1 , (18) 

are independent of the properties of the transition states; 
then, via periodicity, one concludes that 

3=0 

In summary, one might be tempted to regard the load- 
dependence expressions Ijl5|l and H16|l merely as phe- 
nomenological expansions in powers of Fd/ksT; how- 
ever, the crucial feature lies in the physical interpretation 
of the load distribution vectors 8^ and the compliance 
matrices ry^ which, in turn, yields the constraints em- 
bodied in and fljfl . 

4. MOTILITY OF KINESIN UNDER 
VECTORIAL LOADING 
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FIG. 5: Contours of the velocity, V(F X ,F Z ; [ATP]), of kinesin 
at saturating ATP in t he (F Xr F z ) plane as derived from the 
data of Block et al. [Tolll Rll2^| . The dotted, dashed and dot- 
dashed loci, related to the microtubule buckling experiment 
of Howard and coworkers [T^l (see Fig. [3J, are discussed in 
the text. 



Our primary aim now, as an application of the forego- 
ing analysis, is to reconsider the GMBH buckling exper- 
iment . As explained in the Introduction with the aid 
of Fig. |3| the crucial feature of that study was the direct 
determination of the perpendicular force F± = F z and, 
independently, the longitudinal or parallel component, 
given by Fn = —F x (since the load is always resistive 
from the perspective of the track). We will suppose, in 
the absence of contrary evidence, that the microtubule 
(MT) does not twist and that the kinesin, when it at- 
taches to the MT, binds on to the nearest protofilament 
which, again we suppose for simplicity, does not spiral 
around the MT Q] . Then any externally generated trans- 
verse force components may be neglected: i.e., F y = 0. 

In order to analyze the buckling data theoretically 
following GMBH, one needs an explicit expression for 
V(F X ,F Z ; [ATP]), the velocity of the motor along the 
track as a function of F x , F z at fixed [ATP] . In default of 
a better description, GMBH postulated a simple linear 
dependence of V on F x and F z . However, we may hope to 
do better by using the recent data of Block and coworkers 
0, 0] (BASL) who imposed assisting and resisting (as 
well as sideways) loads. The BASL experiments stud- 
ied squid kinesin whereas GMBH employed kinesin from 
bovine brains: undoubtedly this and other disparate ex- 
perimental details should play some role. Nevertheless, 
experience suggests that the quantitative effects should 
not be great. 

Accordingly, in Fig. [3] we present two-dimensional 
velocity-contour plots of conventional kinesin in the 
{F X: F Z ) plane at saturating ATP concentration by us- 
ing the simple (N = 2)-state kinetic model as fitted to 
the BASL data. Details of our analysis will be discussed 
elsewhere |l5ll28| . For completeness, however, we explain 
briefly here how these contour plots have been derived. 



In terms of the general N = 2 expression 5] 

V = d(u ui — w wi)/(u + Mi + wo + w%), (20) 

the dependence of the velocity as a function of [ATP] un- 
der (effectively) zero load may be accounted for, following 

ug = *g[ATP], Wq = k' [ATP]/(l + [ATPj/co) 1 / 2 , (21) 

in which the final, square-root factor (introduced to take 
account of an ATP regeneration system and a related 
trend in the stall force) plays only a small role. Then, 
with d = 8.2 nm a good description of the kinesin data 
of BASL is given by the rates 

fcg = l^AiM-V 1 , w^ = 5.0s-\ u^lOOs" 1 , 
k' = 2.04 x 10~ 3 /dVrV 1 , c = 20 /LtM, (22) 

which, apart for the larger value of k' , are quite compa- 
rable to the original fits [g based on the data of Visscher 
and coworkers HH] (restricted to F x < 0). 

Now Block et al. 01 measured the velocity V only 
as a function of the parallel component F x (even though 
for both F x < and F x > 0). However, in order to use 
(|20|l with the .F-dependence expressions (|15H and HIGH , 
one also needs to know the perpendicular component F z . 
The route to unraveling this puzzle depends, as indicated 
in the Introduction, on a consideration of the geometry of 
the bead-tether-motor-track configuration as illustrated 
in Fig. |21 However, further properties of the motor and 
the experimental set-up may need to be accounted for. 
Concretely, we desire a formula, say F z = T Z (F X ), that 
correlates the imposed (but unmeasured) perpendicular 
component F z with the observed parallel component F x 

0. 
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The most basic hypothesis — termed Model 15, 28] 
— is to suppose, following Fig. [21 that the tether angle 6 
simply switches sign when F x passes through zero: this 
is expressed by 



Model 0: 



F z =F X cot G(F X ) = cu\F x \, (23) 



where cm = |cot0| follows from (j2J). Reasonable values 
for beads of diameter 2R = 0.50 /im, a fluctuating scale 
Az = 5 nm [12j and a tether length / = 60 nm yield 
8 ~ 35° and c\\ ~ 1.44. 

In practice, the thermal fluctuations already recog- 
nized by allowing for Az in Fig. [21 come into play strongly 
when \F X \ < 0.5 pN. Accordingly, a more realistic hy- 
pothesis (which is tested further in Sec. below) is em- 
bodied in the smoothed form 



Model I: 



F z = T Z (F X ) = c n / F* 



(24) 




in which the fluctuation amplitude Fq — 0.3 pN proves 
appropriate in light of observed force fluctuations |3, El • 
Study of the BASL data, however, reveals an unexpected 
and strong asymmetry in the relation F z — F Z {F X ) for 
kinesin 15. 281. This can be represented effectively via 
an additive term as 



Model ILK = 



The amplitude Fx, which is found to be close to 2.0 pN, 
measures the asymmetry. Although the new term may 
appear as an artificial construct it can be interpreted 
rather directly in mechanical terms 0,H3- Note, in- 
cidentally, that the factor in large parentheses merely 
represents a smoothed Heaviside step function vanishing 
rapidly for F x < 0. 

On this basis, successful fits to the BASL data [in- 
cluding measurements of the randomness r(F x ; [ATP]) = 
D/dV] are achieved with the load distribution vectors 

6+ = (0.120, 0, -0.043), 8+ = (0.032, 0, -0.026), 
8q = (0.950, 0, 0.100), 0f = (-0.102, 0, -0.031). 

(26) 

Of course, these satisfy (|17f) so that only six independent 
fitting parameters are entailed. The compliance matrices, 
t/q and t]q , do not need to be invoked although one might 
reasonably presume that they are not totally negligible 
in reality. 

Now we may note, using © and (0J, that these load 
distribution vectors imply a very small forward step of 
only do — 0.1 to 0.2 nm (allowing for the fitting uncer- 
tainties) on binding ATP but a much larger downwards 
displacement, namely Azq = —0.5 to —0.8 nm. In this 
sense, then, the motor appears to "crouch" before it com- 
pletes the main forward step of magnitude d\ ~ 8.0 to 8.1 
nm. It should also be remarked |g, 0] ■> that the relatively 
large positively directed value of 9q means that it is the 



reverse rate wo that changes most under a resisting load 
and thereby leads to the motor stalling (V = 0). 

With the parameters presented above, the contours 
plotted in Fig. follow from (121 )|) . We may remark that 
when the level of ATP is lowered the form of the con- 
tours changes little qualitatively but the scale is reduced 
in accord with (|20|I - H22|I which imply close-to Michaelis- 
Menton variation of V([ATP]) at zero load. It should also 
be noted that the direct sampling of the (F„, F x ) plane by 
the Block and coworkers experiments |3l S El 03 * s essen ~ 
tially confined to the locus specified by the F z = F Z (F X ) 
relation and so, by Model I, is roughly given by F z = 
c \\\Fx\ > 0-3 pN. In principle, the behavior further from 
this locus could be somewhat different than portrayed. 
In particular, the fact that the fitting was achieved with- 
out explicitly employing any compliance coefficients is 
responsible for the steep rise of the velocity contours in 
Fig. [SJas the stall force, Fs ~ — 5.9 x pN, is approached. 
Thus, in the absence of the compliance matrices and fur- 
ther ^-dependent corrections, the stall condition V = 
combined with the periodicity constraint l|17|) for the load 
distribution vectors implies independence of F y and F z . 
Some evidence on the likely magnitudes of the compli- 
ance effects is available from analysis of the transverse 
loading experiments of Block et al. [Tot l2£jj . Shifts in the 
contours, away from the T Z {F X ) locus, as large at 10 to 
25% might be realized. 



5. BUCKLING OF A MICROTUBULE BY A 
MOTOR 

As explained above, in the experiments of GMBH 0] 
the minus end of a microtubule was clamped to a flat 
glass substrate on which kinesin molecules were sparsely 
tethered. When an individual motor encountered and 
bound to the thermally flexing MT, it moved towards 
the plus end and soon buckled the MT: see Fig. |3J Then 
in time sequence (at 30 frames per second), the succes- 
sive shapes of the MT were recorded. From the progress 
of the motor along the contour of the flexing MT the 
time varying velocity, V(t), could be found by fitting 
the observed displacement curves. This was plotted (see 
|13|) against the corresponding force components, F x and 
F z , derived from the shapes. Thereby, GMBH gener- 
ated displays somewhat resembling those shown in Fig. 
|5J although with significant fluctuations (and digitizing 
noise) in all three variables V, F x and F z . In the cases 
reported, the parallel components, = — F x , started 
close to the stall force Fs ~ 6.5 pN and decreased more 
or less steadily with time after buckling, while the per- 
pendicular components, Fx = F z , increased from close 
to zero. In parallel, the velocities started from low val- 
ues, close to stall, but increased, sometimes very steeply 
with changing force, and eventually, for F± and Fu in 
the range 2—5 pN, passed through maxima as high as 
850 to 1100 nm/s. Since these velocities correspond or 
exceed the highest normally seen in standard single-bead 
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FIG. 6: Velocity-force plots mimicking those observed in the 
microtubule buckling experiments of Howard and coworkers 
fl^l . These plots are based on the putative force trajectories, 
A and B, respectively, shown in Fig. |3 The open circles are 
for the perpendicular force components, F± = F z , the solid 
circles for the parallel components, Fu = —F x . 



assays under low loads |-F X | < 1 pN, GMBH concluded 
that the primary action of a perpendicular load was to 
increase V. 

The (V, F XiZ ) plots in Fig. [SI which roughly mimic the 
GMBH observations, were generated from the veloc- 
ity contours shown in Fig. [3] by postulating the dashed 
and dot-dashed force trajectories labelled A and B, re- 
spectively. While our (V,F XjZ ) plots resemble those of 
GMBH, they attain much lower maximal speeds. This 
is clearly a consequence of the facts (a) that, as normal, 
V(F X = F z = 0) ~ 800 nm/s and (b) that for V > 200 
nm/s the contours of V in Fig. [S] slope upwards to the 
right so that, contrary to the conclusion of GMBH, the 
predominant effect of increasing F z is to reduce the ve- 
locity. It should also be remarked that the close to 45° 
slope of the velocity contours for F x > directly reflects 
the observations 0, 0] that even large assisting loads 
do not result in any significant increase in V. 

Before discussing possible reasons for the disagreement 
with the buckling experiments, we ask what form the cor- 
responding trajectories should take in the (F x , F z ) plane. 
We will learn that the putative trajectories A and B in 
Fig. [S] are not very plausible. 

Now, the persistence length of a microtubule is about 
6 mm Thus in analyzing their buckling data |l3j| . 



GMBH considered an MT as a stiff, uniform rod with 
a measured flexual rigidity EI (where E is the Young's 
modulus and I is the moment of inertia of the cross- 
section). The theory of bending a uniform rod is well 
established |29j; but we will sketch it briefly in order 
to understand how the force trajectory in the buckling 
experiment may be derived. This then enables one to 
calculate the velocity of the kinesin molecule as it moves 
along the MT under the resulting induced load. 

We will now assume that the MT is clamped at the ori- 
gin of the (x, z) plane, while the motor is located close to 
the x-axis, say at (Lq, 0), where Lq is the initial distance 
along the MT to the binding site in the clamped position: 
see Fig. |21 Note that the coordinates x and z here dif- 
fer from those introduced in Fig. ^ where the x-axis was 
taken parallel to the MT axis. Similarly, we will always 
denote the parallel and perpendicular components of the 
force relative to the MT as F\\ and F±, respectively. 

When a two-component force, F = (Fu , Fx), is exerted 
on the MT by the tethered kinesin at (Lq, 0), the shape 
of the buckled MT (considered as a uniform rod) satisfies 
the beam equation which can be written 



g = -/3 2 sin[^) 



tp F ] with (i 2 = F/EI, (27) 



where s is the arc length measured along the MT from the 
clamped origin (0, 0) while 0(s) is the tangential angle of 
the rod at the point [x(s), z(s)] with respect to the x axis; 
in addition, ifp denotes the angle of the applied force F 
(with respect to the x axis) at the location of the kinesin 

which defines the pivot point, while F = y' ' F 2 + ~F~f_ is 

the magnitude of the buckling force. 

Assuming that the clamping of the MT is tight and 
that no torque is exerted at the pivot point, the solution 
of this equation must satisfy the boundary conditions 



(s = 0) = 



(de/ds) s 



0. 



(28) 



where L is the total (time varying) arc length to the pivot 
point from the origin at the clamped position. Further- 
more, the pivot point is fixed at the initial location of 
kinesin, which itself does not move spatially, so that one 
must have 



x(s = L) 



/ cos 9(s)ds = Lq, 
Jo 



z(s=L) = 



sme 



= 0. 



(29) 



Now, if the angle 6(s) remains small — as it will when 
the MT starts to buckle — one may accept the approx- 
imations (dx/ds) = cos 9 w 1 and (dz/ds) — sinO w 6 
which lead to 

6 « (dz/dx), (d 2 6/ds 2 ) w (d 3 z/ds 3 ) « (d 3 z/dx 3 ). 

(30) 

Expanding (|27|l up to first order then yields the linear 
equation 



d 3 z ^ 02 dz 
dx 3 dx 



(31) 
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for the displacement z{x) of the MT from the x axis. 
This is subject to the boundary conditions z(0) — z'(0) = 
z"(0) = z(L ) = which then give 

z(x) = <pf(x — Lq + LqCOs(3x — (3~ 1 sinflx). (32) 

But by (|29|l this form implies that the buckling force, 
given by F — (3 2 EI, must satisfy the equation 

tan (3L = (3L ~ 4.493. (33) 

Note that in this small-0 approximation the magnitude 
of the force remains constant at F* ~ EI(4A93/ Lq) 2 as 
the rod initially buckles. This in turn leads to a force 
trajectory in the (F X ,F Z ) plane that is just a semicircle 
of radius F*. Of course, the small-angle approximation 
fails when Fx increases. 

Exact solutions to the nonlinear equation ((27(1 can be 
expressed via an elliptic integral Q, El as 



J<Po v 1 — At sin (j> 

in which the modulus k and elliptic angle </>(s) are related 
by 

ksm.(j>{s) = sin \[9{s) — tpp], k 2 = sin 2 ~ <Pf)i 

(35) 

while 9l = 6(s — L). Note that the solution represented 
by 1)34(1 already satisfies the boundary conditions ((28(1 . 
However, the initial elliptic angle 0o and the modulus k 
must be found so as to reproduce the correct pivot rela- 
tions 1(29(1 while /3 oc y/~F determines the scale of the arc 
length. The parallel and perpendicular force components 
for an MT buckled to a contour length L > Lq are then 
given by 

F\\ = Fcos(6 L - tp F ) = -F(l - 2k 2 ), (36) 
Fx = Fsm{6 L - <Pf) = 2Fk\J\ - k 2 . (37) 

Finally, for a given L > Lq, one may integrate 1(29(1 
numerically using ((34|) and l(35|l to obtain the actual tra- 
jectory of MT buckling in the (F x , F z ) plane. The dotted 
curve in Fig. [5] represents such a trajectory when the ini- 
tial parallel force is set equal to the stall force F$ ~ 5.9 
pN. For small Fx = F z , the trajectory approaches the 
circle described by the small-angle approximation. 

The (V,F X _ Z ) buckling plot following from the calcu- 
lated force trajectory (dotted curve in Fig.[5J is presented 
in Fig. Evidently, this differs significantly from the re- 
sults of GMBH and from the postulated forms shown 
in Fig. Indeed, the key experimental finding, namely, 
that the velocity, although noisy, exhibits a maximum 
above 850 nm/s at intermediate force levels is not re- 
produced: on the contrary one finds V < 650 nm/s in 
this region. On the other hand, the theory presented by 
GMBH, which simplifies the functional form of the ve- 
locity to a linear expression in the force (see 01 Fig- 9c) 
generates qualitatively similar behavior to Fig. \7\ Thus 



800 




2 4 6 



F x , z ( P N) 

FIG. 7: A velocity-force plot for a microtubule buckling ex- 
periment, as in Fig. |S| but based on a force trajectory calcu- 
lated using theory for the bending of a uniform slender rod 
and shown as a dotted curve in Fig. 

Fh falls steadily while V always increases and Fx passes 
through a maximum: see 0] Fig. 10. 

In calculating the force trajectory for buckling an MT, 
we assumed that the clamping is tight and that no torque 
is applied by the motor protein at the pivot point. These 
assumptions seem reasonable based on the fits to the ex- 
perimental data presented in Fig. 5. Nevertheless, as 
one can detect in this figure, the experimental data, espe- 
cially at the beginning of the buckling event, show small 
deviations from the fitted curve near the pivot point; in- 
deed, the microtubule appears to be slightly bent at the 
pivot point. This suggests that the motor protein may 
actually exert a torque on the MT which might lead to a 
rather different force trajectory. 

Other possible sources of uncertainty in the experi- 
ments and their interpretation were discussed critically 
by GMBH, including clamp looseness, misalignment of 
the motor and the axis of clamping (z = in Fig. |3J), 
protofilamcnt number variations and slight, preformed 
bends in the microtubules, etc. Nevertheless, it is diffi- 
cult to understand how the large maximal velocities ob- 
served, exceeding 900 nm/s, could be significantly in er- 
ror. If this conclusion stands — and, clearly, experiments 
in which F z and F x can be controlled more directly and 
with fewer uncertainties are to be desired — a satisfac- 
tory explanation remains to be found. It is certain that 
the interaction of the heads, i.e., the motor domains with 
the microtubule plays a crucial role in kinesin motility: 
indeed, a recent normal mode analysis of simple protein 
network models j20] highlights this feature as a distinc- 
tion separating kinesin from standard myosin and the 
Fl ATPase motor. It seems possible, therefore, that the 
stressed state of the microtubule in the buckling experi- 
ments could be a prime cause of velocity enhancement. 
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TETHER LINKAGE UNDER LOAD 



Underlying our analysis |14l Il5l |28| of the Block et al. 
experiments imposing reverse and transverse loads on ki- 
nesin 0, are the expressions embodied in Models I 
and II for the perpendicular force, F z — !F Z (F X ), induced 
via the tether and bead by the longitudinal (or parallel) 
load: see ij^ljl and (|2*5l) in Sec.^] It is clearly of interest to 
attempt to test these forms against experiments that ex- 
amine the transmission of force in the bead-tether-motor- 
track assembly. To that end we may utilize data obtained 
by Svoboda and Block |l2( who studied the elasticity of 
the bead-kinesin-microtubule linkage. 

These authors measured the ratio of the velocity V\> of 
a bead trapped by optical tweezers to the velocity V s of 
the stage to which the MT was attached, in the presence 
of the nonhydrolyzable ATP analog AMP-PNP. This pro- 
duces a rigorlike association between the kinesin motor 
and the MT. The ratio may then be expressed in terms 
of elasticities as 



V h /V s = K mot /(K mot +K h 



(38) 



where K mot and Kb are the stiffness of the bead-motor- 
track linkage and of the bead in the optical trap, re- 
spectively. The ratio was measured as a function of the 
displacement Xb of the bead from the center of the trap. 
The individual observations at a laser power of 62.5 mW 
are presented in Fig. El Although the data are noisy ow- 
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FIG. 8: The ratio of bead velocity, Vb, to stage velocity, V s , 
as a function of bead displacement, x^, from the trap center, 
as observed obtained by Svoboda and Block fl^ . The solid 
circles are averages of the raw data (pluses) using a 10 nm 
bin size. The dotted lines are fits of Model I for the relation 
F z — J- Z (F X ), to the solid circles with noise levels Fq = 0.3 
pN and 0.5 pN, while the solid lines are corresponding fits for 
Model II. 

ing to the Brownian motion of the bead and to linkage 
heterogeneity [lj] , the mean ratio rises slowly as the dis- 
placement increases from the smallest observed value at 



Xb = 50 nm: see the solid circles which have been cal- 
culated by averaging the data binned in intervals of size 
10 nm. We may then ask how well our models fit these 
data. 

To proceed, note that the longitudinal force, F x , ap- 
plied by the optical trap for a bead displacement Xb sat- 
isfies 



-F x = K h x h = K mot Al x , 



(39) 



where AZa; is the change under tension of the projection of 
the kinesin tether on to the MT or x axis. By assuming, 
as is reasonable lj, that the kinesin tether is rigid, we 
have Al x — I sin 1 9 1 where I is the length of the tether 
while 8, here negative, is the angle of the tether from the 
vertical axis: see Fig. |3 After some algebra, one finds 



Vs 



x b 



Xb 



' sin 161 ' 



sin O = 



F x 



(40) 



Finally, we may use (|39|l in combination with the mod- 
els (I24|l and l|25|l for J- Z (F X ), to predict the ratio as a 
function of Xb- For the data in Fig. ODone has Kb = 0.03 
pN/nm H we take I = 40 nm, the models I and II 
yield perfectly satisfactory (if not very informative) fits 
to the data as evident from the solid and dashed lines 
in the figure. Note the two values assumed for the force 
fluctuation F in and (|25J| • The tether length fitted 
here is shorter than what would be expected for a free ki- 
nesin molecule on the basis of photomicrographs , say, 
I < 60 nm. However, when a bead is chemically bound 
to a kinesin molecule it seems likely that some length 
of the tether may also be attached to the surface of the 
bead. Thus there are no grounds here for questioning the 
adequacy of the tether-based models for J r z (F x ). 



7. SUMMARY 

We have extended the previous simple mechanochem- 
ical kinetic models for motor protein motion 0, El t° 
accommodate a three-dimensional free energy landscape 
for the point of attachment of the tether to the body of 
a motor that moves under a vector load F. The load- 
dependence of the various forward and reverse rates de- 
scribing the mechanochemical progress of the motor as it 
takes a single overall forward step of size d, are embodied, 
to first order in F, in a set of load distribution vectors 
0+ and 6j . These, in turn, relate directly to the forward 
substeps, dj(< d), taken by the motor along its track 
and also serve to localize the corresponding intermediate 
transition states. In quadratic order in the components 
(F x ,F y ,F z ) of F, a set of compliance matrices, rf^ and 

rjj , comes into play. 

Even when perpendicular (or vertical) force compo- 
nents, F z , are not purposefully imposed on a motor, as 
in standard single-bead optical-trap assays, consideration 
of the geometry of the bead-tether-motor-track configu- 
ration demonstrates that loads with F z > are induced. 
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Furthermore, in switching between resistive and assisting 
loads (parallel to the track) accounting for these vertical 
components proves imperative. 

The general analysis has been used to reconsider 
the microtubule buckling experiments of Howard and 
coworkers 13] in which the values of F z were measured. 
To that end the approach was first applied to the recent 
experiments of Block et al. 10] in which, in particular, 
assisting loads imposed on kinesin were found to have 
little if any accelerating effect. The resulting fits provide 
velocity contours in the (F X ,F Z ) plane for kinesin under 
specified [ATP] : these are needed to address the buckling 
data of Howard and collaborators. Incidentally, however, 
the analysis also indicates that, on initially binding ATP, 
a kinesin motor "crouches," i.e. moves downwards closer 
to the microtubule, prior to hydrolyzing ATP and step- 
ping directly forward by ~ 8 nm. Further and more de- 
tailed discussion of the Block et al. exp eriments on ki- 
nesin [l(| is provided elsewhere 15, 28]; but, as shown 
here, the associated modelling of the transmission of ten- 
sion via the bead-tether-kinesin-microtubule linkage is 
consistent with previous measurement by Svoboda and 



Block 12]. 

On the basis of the buckling experiments ^3 Howard 
and coworkers concluded that a perpendicular force F z = 
3-5 pN, together with a longitudinal resisting force, \F X \, 
of similar magnitude, results in kinesin velocities in ex- 
cess of 900 nm/s. Our analysis does not substantiate this 
inference suggesting instead that velocities no larger than 
say 600 nm/s or so, should have been observed. A reso- 
lution of this challenging discrepancy might reside in the 
influence of microtubule stress or curvature on kinesin 
motility. 
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